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ABSTRACT 

“Astronomical” or “circumstellar” silicate optical functions (real and imaginary indices of 
refraction n(A) and fc(A)) have been previously derived from compositionally and structurally 
disparate samples; past values were compiled from different sources in the literature, and are 
essentially kluges of observational, laboratory, and extrapolated or interpolated values. These 
synthetic optical functions were created because astronomers lack the quantitative data on amor¬ 
phous silicates at all wavelengths needed for radiative transfer modeling. This paper provides 
optical functions that (1) are created with a consistent methodology, (2) use the same sample 
across all wavelengths, and (3) minimize interpolation and extrapolation wherever possible. We 
present electronic data tables of optical functions derived from mid-ultraviolet to far-infrared lab¬ 
oratory transmission spectra for two materials: iron-free glass with chondritic/solar atmospheric 
abundances, and metallic iron. We compare these optical functions to other popular n, k data 
used to model amorphous silicates (e.g., “astronomical” or “circumstellar” silicate), both directly 
and in application to a simple system: the dust shell of the post-AGB star HD 161796. Using 
the new optical functions, we find that the far-IR profile of model SEDs are significantly affected 
by the ratio of glass to iron. Our case study on HDi61796 shows that modeling with our new 
optical functions, the mineralogy is markedly different from that derived using synthetic optical 
functions and suggests a new scenario of crystalline silicate formation. 

Subject headings: (ISM;) dust, extinction - methods: laboratory - ultraviolet: ISM - infrared: ISM 


1. Introduction 

Silicate grains dominate the emission and ex¬ 
tinction processes in many astrophysical envi¬ 
ronments but astronomers still do not know 
exactly what composition or structure consti¬ 


tutes these dust grains or how these proper¬ 
ties change from location to location. As- 
tromineralogy, the study of the precise and 
detailed nature of dust grains in space, has 
developed rapid ly over the past few decades 
(see reviews in ISneck et al.l 119971 : ISneckl Il998 : 
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Sueck et al.ll200(Tl: MolsteJ200r)l: Waters & Molster 

1999; 

Hennind 2003; Kwok 2004; Pitman et al. 

2010; 

Gnha Nivosri et al. 2011al: Sueckl I2OI2I). 


However, interpretation of dust grains is still lim¬ 
ited by the paucity and quality of mineral data 
used for comparison with and modeling of astro¬ 
nomical observations. 


While dust in space is particulate in nature, 
there is some discussion about whether it is ap¬ 
propriate to use bulk materials for spectroscopy of 
cand i date solids fe.g. ISoeck. Hofmeister fc Barlow 


19991: E ennindI 2 OOIIL and references therein). It 
has already been shown that thin films (slabs) 
and dilute particle di stributions result in the same 
spect ral features ( Speck. Hofmeister fc Barl^ 
I 999 II . Moreover, having optical constants from 
spectra of bulk materials is important because 
they set an upp er limit for extinction t hrough a 
certain material ( Hofmeister et al. 200911 . 


Although laboratory measurements of optical 
functions for astronomical applications are of high 
value, they characterize individual analogue ma¬ 
terials and are not the real dust we see in space. 
However, good quality extinction properties of the 
components of the real space dust can be derived 
from these analogue materials. Likewise, deriving 
a set of effective optical functions from astronom¬ 
ical data would be acceptable if the temperature 
and density structure, and grain size and shape 
distribution of a system was well known. Here we 
aim to provide optical functions that allow us to 
determine these details in astronomical environ¬ 
ments more completely. 


There are two basic approaches to identifying 
dust minerals in space: (1) Work backward from 
astronomical observations where temperature and 
density structure of a system is well known to de¬ 
rive effective optical functions that approximate 
those of real dust in space. (2) Work forward from 
laboratory measurements of dust analogue min¬ 
erals whose compositions and impurities are well 
characterized to derive optical functions of miner¬ 
als predicted to occur in space. 


Ideally, we should use both approaches and find 
where (1) and (2) match, however, astronomers 
do not take the dual approach. Approach (1) 
alone does not reveal anything about the chem¬ 
istry; those optical functions are not useful for 
diagnosing changes from location to location and 
should not be used to interpret compositional ef¬ 


fects elsewhere. We are taking approach (2) by 
providing dust analogue data that will provide in¬ 
sight into the chemistry. Moreover we are correct¬ 
ing past efforts which made kluges that hybridize 
rather than compare the two approaches. 

The simplest approach for determining the na¬ 
ture of dust grains in space is to match the posi¬ 
tions, widths, and strengths of observed spectral 
features with those seen in laboratory spectra. 
This is usually achieved by fitting a continuum to 
the observed spectrum and dividing (or subtract¬ 
ing) to derive emissivities of the observed spectral 
features. These emissivities are then compared to 
absorption spectra0 produced in the laboratory. 
Sometimes the emission/absorption profile is cal¬ 
culated from optical functions (complex refractive 
indice s: n\ and kx; see e.g. lOnha Nivogi et al. 
2011al ld. Radiative transfer (RT) modeling uses 
the optical function^ minerals to predict how a 
given object shoul d look both spectroscopically 


and 

in 

imaees ('e.e. IPinte et al.l 

2008; 

Eisner et al. 

2005 


Meixner et al. 

2002; Men’shchikov et al. 

2002 

; Voors et al.l 2000 

: Sorrell 1990l) 

RT model- 


ing can be used to determine the effects of grain 
size distributions and mineralogy on the expected 
spectrum, and places constraints on the relative 
abundances of different grain types in a dust shell. 
However, optical functions used in such RT mod¬ 
els must cover a broad wavelength range, from the 
ultraviolet (UV) to the infrared (IR). Deriving op¬ 
tical functions over this broad range requires lab¬ 
oratory transmission or reflectance spectra from 
multiple detectors; this is one of many reasons 
that n, k are hard to find. As a result, various 
compilations have been produced that kluge to- 


^ Typical absorption spectra are derived from transmission 
experiments which yield the absorbance a or the absorp¬ 
tivity A. The details of these dif ferent parame ters and how 
they are derived are reviewed in ISpeckI l(203)' 

^The complex refractive indices of a substance are often re¬ 
ferred to as the “optical constants.” Because these values 
are not constant, but rather vary with wavelength, temper¬ 
ature, and polytype (structure) as well as chemical compo¬ 
sition of the sample, we use the term “optical functions” in 
this work. 

®while RT models use absorption and scattering, or ex¬ 
tinction efficiencies of candidate dust grains, many are 
equipped to use optical functions as inputs and do the 
conversion themselves, then generate the required cross¬ 
sections/efficiencies. Providing optical functions allows for 
more versatility in the modelling (e.g. grain size and shape 
effects) 
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gether observed and laboratory-collected data in 
order to produce optical functions to fill the full 
wavelength grid required by RT codes. 


Still the most widely u sed optical fu n ctions 
for silicate s are those from l^aine fc Led (1984 ) 
(IDLI 1984, hereafter), Draine (llQSd 12003 ) and 
lOssenkopf, Henning, fc Mathis ()l992ll (foHMlll9^ , 
hereafter) which were originally generated 30 years 
ago and included kluges out of necessity. How¬ 
ever, we are now able to provide optical func¬ 
tions which are kluge-free from the far-IR to the 
UV (0.19 < A < 250 ^m). Here we present new 
laboratory-generated complex refractive indices 
for an iron-free glassy silicate of chondritic/solar 
atmospheric composition using a single sample 
for all wavelength regimes and compare it to the 
most commonly used “astronomical silicate” op¬ 
tical functions. We apply our new silicate glass 
optical functions to modeling the dust around 
HD 161796, an oxygen-rich pre-planetary nebula, 
in order to compare our new optical functions 
to those most commonly used. In addition we 
present new optical functions for metallic iron in 
order to properly model the relevant astrophysical 
environments. 


There are a number of published lab¬ 
oratory datasets on metallic iron, (see also 
http://www.astro.spbu.ru/JPDOC/f-dbase.html), 
however, most have a very limited wavelength 
range. In astronomy, t he mo st widely used data 
are from lOrdal et al.l (|l988h which have both 
low resolution and do not extend to wavelengths 
shortward of t he red -end of the visible spectrum. 
Pollack et al. ( 1994t l also present a figure which 


contains a broad wavelength range optical function 
dataset for metallic iron, but this is compiled from 
multiple samples and thus not discussed further 
here. Likewise, there i s a set of optical functi ons 
for m e tallic iron from iLvnch fc Hunted (jl99ll ) in 
PalikI ( I99l[l . but it also is derived from a com¬ 


pilation of samples. Iron optical properties are 
also important for remote sensing in planetary 
science, but those studi es are usually lim ited to 
the NUV-NIR range Isee lCahill et al.lf20I2L whose 
own excellent data only covers 0.16-1.7/rm). 

The new Fe optical functions improve upon 
the previous datasets because (a) most previous 
datasets do not extend beyond the red-end of the 
visible and (b) many are compilations from dis¬ 
parate samples. This leads to artifacts, such as a 


discontinuity in R at 600 nm. In contrast, our new 
data has the higher spectral resolution, broader 
wavelength coverage, and is derived from the same 
sample through that wavelength range. 


2. Past work: “Astronomical” and “cir- 
cumstellar” silicates 


Several different synthetic optical functions 
were derived at a time when there was insufficient 
laboratory data to be able to model astrophysical 
environments across a broad wavelength regime. 
A detailed discussion of lab data for Si-based ma¬ 


terials can be found in IColangeli et al.l (|2003ri as 
well as ISneck et al.l pOlll l. and a wide r discus¬ 
sion o f cosmic silicates ca n be found in IHenning 
( 20101) . Speck et al. ( 201ll) made detailed compar¬ 
isons between previous published laboratory data 
on silicate glasses. 

Compared to the wealth of data on poten¬ 
tial astrominerals, our literature search (e.g. 
http://www.astro.sp bu.ru/JPDOC/f-dbase.html 
Henning et al.l Il999l and references therein) re¬ 
vealed that although silicates have been studied 
in more detail than other compounds, there are 
still shortcomings in the available data with re¬ 
spect to silicate glasses. Specifically, existing lab 
data on glasses do not allow astronomical spec¬ 
tra to be interpreted reliably in terms of either 
M^+/Si or Mg/Fe ratio, which are important tools 
for discriminating between competin g dust forma- 
tion models f c.f. ISoeck et al. 2011 . Figs. 7-9). 
Pitman et al.l ( 2013lj) compared the UV-Visible 
extinction coefficients (and imaginary index of 
refraction; k) for a range of silicates, both amor¬ 
phous and crystalline with those from a variety 
of previously publish ed samples. The figure from 
Pitman et al.l ( 2013bh (reproduced in Fig.[T]) illus¬ 
trates that k values vary wildly for both crystalline 
and amorphous silicates in the UV-vis, whether 
the datasets are from the lab or from model- 
ers. The 70% and 80% Mg-rich silicate glasses of 
Dorschner et al.l (|1995l) . las er-ablated amo r phous 
forsterite and enstatite from IScott fc DulevI (119961) 
and the “as trono mical silicates” of , ^ _( 1984[ ). 
OHM ( 1992 ). and David fc Pegourie ( 1995 ) have 
a similar u pper limit on k. The k values for the 
remaining iDorschner et al.l ( 1^95 ) Mg-rich sil i- 
cate glasses, sample 1-S from J^ager et H] (1 1994 ) 
and enstatite and olivine from lEgan fc Hilgeman 
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(Il975a ) a re a fa ctor of ten less. The datasets from 


HufFmanI (1975) are in a class by themselves (set¬ 


ting the uppermost limit on A:), as are the Fe-rich 
silicate glasses. It is clear from that ostensibly 
similar compositions give rise to very different op¬ 
tical functions. In addition for crystallinity and 
composition differences, the variation in derived 
optical functions is also due to the use of powders, 
rather than bulk samples. Studis involving pow¬ 
der samples often overestimate the k values be¬ 
cause of internal reflections and scattering. This 
problem also afflic ts bulk samples which i nclude 
impe r fections (see Hofmeister et al] 20091 : ISoeck 
2014 ; ISoeck. Hofmeister fc Barlowlll999 , for more 


detailed discussion of this issue.). 



Fig. 1.— NUV-NIR Imaginary index of re¬ 
fraction, k, (a.k.a. extinction) of commonly 
used “Astronomical ” silicates (reproduced from 
Pitman et al. 2013b[) 


Here we focus our discussion on the most 
commonly used datasets for analysi s of s ilicate 


my _ _ 

dust dPH Il984ll . iDrainel (Il985l 1200311 and lOHM 


(jl992n . Other artificial “circumstellar,” or “astro¬ 
nomical” silicate s, based in part on real mineral 
data, also exist 

David fc PegQurie 



1991: Laor fc Draine 1993: 


199511 . The optical func tions in 


Suhl (|199lll were adopted from |DL| (1198411 optical 
functions short ward of A =8 um and are thus not 
discussed here. IPavid fc Pegourie ( 1995ll followed 
a similar methodology to lPlJ ( 1984[ l. but using the 
then new IRAS observations of ~300 red super¬ 
giants. However, it is not widely used, prob ably 
becau se lOHMI (1199211 had alrea dy supp l anted IdlI 
( 1984 1 for some purposes, and lOHM ( 1992 1 in¬ 


corpo rated t he source data forlDavid fc Pegourie 
( 199, 5l l (from Davi^^JPanoular 1990t l into their 
dis cuss i on. Ibaor fc Drainel (|l993 1 is mostly based 


on 


, but with tweaks at the high energy 
end (> 20eV) which are calculated for an assumed 
density and MgFeSi 04 compositi on, but these 
modifications were supplanted by Draine ( 2003| l 
and n ot considered further here. Kim fc Martin 
( 1995 1 found a problem using IDl] ( 1984h to model 
UV polarization which they attributed to the use 
of crystalline olivine data in the derivatio n of 
“astro nomical silicate” optical functions in iDL 
( 1984l l. Thus they proposed using measurements 
from volcanic glass in the UV to alleviate the 
problem. This is essentially still a kluge. In l ight 
of the above, in this work ^ we w i ll foc us on [PL 
( 1984 1: Draine ( 2003ll and OHMI ( 1992ll datasets 
in our discussion. 


Artificial “cosmic,” “circumstellar,” “interstel¬ 
lar,” or “astronomical” silicate, optical and dielec¬ 
tric functions are not optimized at all waveleng ths. 
For ex ample, the “astronomical silicate” from Im] 
(Il984ll is optimize d for the mid - IR, b ut not in 
the U V-vis range (IPollack et al. H194ll . [Sorrel] 
( 1990ll found that while lDil ~( 1984 1 “astronomi¬ 
cal silicate” could produce a satisfactory match 
to the 10 fj,m feature in a protostellar object, it 
coul d not simultaneously match the UV region 
(c.f. Kim. Martin fc HendrvI 1994h . 


Synthetic “astronomical” or “circumstellar sil¬ 
icate” optical functions are often favored because 
they have broad wavelength coverage, especially 
at short wavelengths (UV and X-ray). As stated 
earlier, optical functions derived from merged lab¬ 
oratory spectra are collected piecemeal using mul¬ 
tiple detectors; merging must be performed judi¬ 
ciously by the user. Furthermore, laboratory de¬ 
rived n and k values typically do not cover very far 
into the UV and X-ray range because these detec¬ 
tors must be hand built or, if commercially avail¬ 
able, are very expensive. However, the trade-off in 
using synthetic “astronomical” or “circumstellar” 
silicate optical functions is that chemical composi¬ 
tion has not been preserved across all wavelengths. 
The optical functions currently in widespread use 
were constructed by a combination of forward 
calculation from laboratory spectra of chemically 
disparate samples and backward calculation from 
astronomically-observed dust opacities, together 
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with pieces that are theoretically derived as illus¬ 
trated in Figure [2j Those datasets also include 
modifications to specifically match certain astro¬ 
nomical observations. Consequently, whereas the 
synthetic optical functi ons match some interste llar 
spectra (in the case of lDlj^lQ^ Draine 2003 ) or 
circumstellar spectra (for IOH1Vllh992t) , those n and 
k values are not useful for diagnosing changes from 
location to location and should not be used to in ¬ 
terpret comp ositio nal effe c ts. In deed, IDLI ( 1984 ). 


Draine ( 2003 ). and OHM ( 19921) explicitly stated 


such caveats, but their warnings have been largely 
ignored because better alternatives have not been 
available. 



Fig. 2.— Source data for different wavelength 
regimes in previo us “astrono mical silicate” opti¬ 
cal functions from DL ( 19841) . Solid lines are the 
n values, dashed lines are the k values. These n 
and k are merged from a mix of experimental, ob¬ 
servational, and theoretical data or assumptions. 
Main figure has linear axes: x-axis is wavelength 
in /im; y-axis is the value of the real or imaginary 
refractive index. Inset is the same as main figure 
but with log-log axes. 


IdlI (19840 produced synthetic optical func¬ 
tions for both “astronomical” silicate and graphite 
in order to model interstellar dust. Bec ause of the 
paucity of laboratory data at that time, Im] (|i984 
compiled data from multiple sources in order to 
achieve continuous wavelen gth coverage. The ful l 
procedure is described by iKim fc Martini (11994 
and summarized in Fig. [ 2 ] 

The updates bv lDrain3 ( 20031) apply mostly to 
the short wavelength/high energy end of the spec- 


^Electronic tables of the dielectric function (e) and equatio ns 
to convert these to n, k were published in IPrainS |T983). 


trum such that the values from iDLl (|l984r) were 
slightly modified: ( 1 ) a feature was removed at 
= 6.5/tm“^ (0.15/tm); ( 2 ) for the sub-mm 
range (250-1100/t m) the value of was varied 
by ±12% following iLi fc Drainel ( 200ll) : (3) above 
30 eV (A < 40 nm) £2 was estimated from atomic 
photo-absorption cross-sections, except near ab¬ 
sorption edges. These high energy dielectric func¬ 
tion values were m erged with original functions 
from Drainel ( 19851) between 18 and 30 eV (40< 
A <70nm). 

ohm ] (I1992 ) foll owed a s i milar method to that 
of IdlI |l984h and Draine ( 2003 ) but with dif¬ 
ferent assumptions, c onstra i nts an d observational 
data. In particular, IOHMI (|l992r) used an aver¬ 
age o£jmdTR^_observations from IRAS com piled 
by Volk fc KwokI (jl988l) . As in IDLI ( I984I) and 


Draini 


_ , it was assumed that grains have an 

MRN grain-size distributior0. They also assumed 
a continuous distribution o f ellipsoid sh apes. For 
the NIR-MIR merge region, IOHMl(ll992h used lab¬ 
oratory data from magnetite (Fe 3 04 ) and metallic 
iron to get high enough opacities. 


Both “astronomical” and “circumstellar” sili¬ 
cates will match some spectral features, and can be 
used for comparing optical depths between differ¬ 
ent dusty environments. However, because those 
n and k values are derived from a combination of 
laboratory and observed data as described above, 
they do not represent real solids across all wave¬ 
lengths and thus cannot be used to determine 
the physical and chemical properties of dust in 
space, or how and why dust varies spatially or 
temporally. Moreover, because in the majority of 
systems, most light is absorbed and scattered in 
the NUV-NIR region, calculations of dust abun¬ 
dances are based primarily on optical properties 
for species that are expected to be rare ( e.g., 
crystalline olivin^ or are approximations (e.g., 
alumina/corundum AI 2 O 3 ). In addition, because 
m etallic iron i s invoked either directly (in the case 
of lOHMl[l994 or indirectl y (using “dirty silicates” 
in the case of lDrainel[2003) , the effect of iron abun¬ 
dances on dust opacities cannot be determined 
with certainty. Because iron is predominantly gen- 


®i.e., the number of grains of size a is given by n(a) propor¬ 
tional to a~'>, where n is the number of the grains in the 
size interval (a, a + da) and q = 3.5 ;amin = 0.005/rm; and 
Umax = 0.25/rm llMathis et. allllF^ . 

® based on the interpretation of observations 
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erated in Type la supernovae and the a-rich silicon 
and magnesium atoms are predomi nantly made 


in co r e-collapse (Type II~) supernova (IBulbul et ah. 
2012 : Iwamoto et al. 19991 : iNomoto et alJ I 2 OO 6 II . 


the relative abundances of these dust-forming el¬ 
ements may not scale together as we move from 
place to place, or especially back in time. The sep¬ 
aration of iron is essential to understanding chem¬ 
ical evolution. 

Finally, while optical functions do not de¬ 
pend on grain size/shape; the resulting absorption 
and scattering properties do. Consequently, the 
derivation of optical functions from (absorption 
and scattering) observations necessarily requires 
assumptions about the grains, in particular their 
size and shape distributions. For all the synthetic 
optical functions above, the grains are assumed 
to h ave M RN size distribut ions. Furthermore, 
IDLI (1984) and Draine ( 2003h assumed a specific 
form of oblate grains, while OHM ( 1992h assumed 
a broad continuous distribution of ellipsoids (in¬ 
cluding spheres, pancakes and needles). 

3. Alternative: “cosmic silicate” glass 

Here we present new optical functions for a 
“cosmic silicate” glass with chondritic/solar at¬ 
mospheric abundances of the major elements ex¬ 
cluding iron. These new optical functions were 
derived from laboratory spectra using the same 
sample of silicate glass from A=0.19 to 250/im. 
T he composition is gi ven in Table [TJ As described 
in Speck et al. ( 201 ill . we excluded iron from the 


sample for two reasons: (1) to better preserve ox¬ 
idation state; and (2) because several theoretical 
models of dust condensation suggest that iron does 
not commonly form silicates, but rather that it 
forms metallic iron grains fe.g.. iLodders fc Feelev 
I 999 II . Iron atoms are almost as common as mag¬ 
nesium and silicon atoms in the solar neighbor¬ 
hood, leading to some debate about the inclusion 
of iron in silicates. However, compositions that 
include grains that are agglomerations of iron- 
free silicates and metallic iron pieces are com¬ 
monly invoked. Indeed, “dirty silicates” were sug- 
gest ed because iron-free si licates are too transpar¬ 
ent ( Merrill fc Stein 1976l l. 

In addition to our n and k for “cosmic silicate” 
glass, we also present n and k for metallic iron in 
order to demonstrate this opacity issue. Because 


our silicate is Fe-free, and we provide Fe optical 
functions, astrophysicists can now confidently in¬ 
vestigate the effect of changing Mg/Fe abundance 
ratios in environments where dust forms. This is 
especially important in lower metallicity environ¬ 
ments, where Mg/Fe should be higher than solar 
even though the absolute number density of atoms 
may be lower. 


Table 1: Composition of “Cosmic Silicate” glass: 

(Nap, iiCap, i2Mgi.s6) (Alp, igSii,85)06 


oxide 

weight percent 

Si02 

54.26 

AI 2 O 3 

4.34 

MgO 

36.58 

CaO 

3.27 

Na 20 

1.61 

total 

100.05 

water 

ppm 

content 


H 2 O 

81 


4. Experimental methods 


4.1. Sample preparation 

4-. 1.1. Glass synthesis 

The sample studied here, “cosmic silicate,” is 
an Fe-free silicate glass with chondritic/solar rel¬ 
ative Mg, Al, Si, Ca, an d Na abundances, syn¬ 
thesiz ed a s described b v fCetson fc WhittingtonI 
( 2007t l and lSoeck et al. ( 201lh bv A. Whittington 
at the University of Missouri in Columbia. The 
melting and quenching techniques used produce 
more consistent, fully amorphized glasses than 
other sample preparation techniques (e.g., chem¬ 
ical vapor deposition, i on-irradiation, lase r abla¬ 
tion; see discussion by ISoeck et al.l 1201 ill . The 
“cosmic silicate” glass is free of bubbles and in- 
tern al crystals; e l ectron microprobe composition 
from Speck etahl ( 2011 1 is given in Table [T] 

Our “cosmic silicate” glass has a trace amount 
(81 ppm) of H 2 O, and shows up as a very weak 
feature around 3 /rm. For a binary, ternary, 
or quaternary silicate glass, we would expect 
to see at least one absorption in the range 
of 2.7 to 4.0/rm, depending on the c o mpos i- 
tion of the glass (jBrown fc Martinsenl Il972tl . 
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Given that e! is proportional to constant +10“^^ 
tim es the concentration at long wavelen gths (see 
e.g. lAndeen. Schnele k. Fontanellalll97'il . we can 
safely neglect the trace H 2 O amount in the IR. 

4 . 1 . 2 . Metallic iron 

High purity electrolytic iron (>99.97%) was 
purchased from Alfa/Aesar@in the form of ~cm- 
sized slabs with ^mm thicknesses. 


4 . 1 . 3 . Sample thickness determination 


We acquired spectral measurements of both the 
“cosmic silicate” glass and the metallic iron sam¬ 
ple from polished slabs (in the UV-VIS-IR) and 
thin films (in the IR only). All polished slab sam¬ 
ples were cut to sub-mm scale with a rock saw, 
then progressively ground, polished with micron¬ 
sized grits by hand, and (for transmission mea¬ 
surements in the IR) were thinned by compres¬ 
sion in a diamond anvil cell to films of below 1 ^m 
thick to achieve well-resolved spectral peaks and 
good signal-to-noise ratios in the spectra. 


Ma ny published papers and iBohren fc Huffman 
( 19831) suggest that grain shape has an effect on 
absorption and scattering cross-sections. Mea¬ 
suring optical properties from slabs, rather than 
loose or dispersed powders, ensures that the op¬ 
tical properties depend only on composition and 
not on grain shape or size. Also, using samples 
with differing thicknesses for different spectral seg¬ 
ments helps to minimize the effect of bac k reflec¬ 
tions ( Hofmeister et al. 20091 : Soeckl 20l3) . 


For the thin films, the thickness of the film is 
the major driver of spectral quality. For the pol¬ 
ished slabs, the quality of the laboratory spectra 
is affected by the quality of the sample’s polish 
(for reflectivity) and the thickness of the sam¬ 
ple (millimeter pathlength for reflectivity, lower 
limit of ~40 /rm achievable in the lab for trans¬ 
mission/absorption). We polished the slab sam¬ 
ples to a surface roughness that is equivalent to 
that of the spectrometer’s reference mirror for the 
IR instrument. In the UV, the high absorbance 
overwhelms the small effect of scattering from the 
surface. 


The thickness of the sample, L, which factors 
into the determination of fc, was measured for the 
polished slabs using a digital micrometer or a cali¬ 
brated reticule in a doubly polished microscope. If 


spectral fringes were observed, we calculated the 
slab thickness from Eq. (U 


L={2nyisAiy) \ ( 1 ) 

where v is frequency in cm“^ and Uyis is a tie- 
point for the real index of refraction in the vis¬ 
ible. Our best determination of Uyis of a thin 
(^mm-scale), doubly polished slab of “cosmic sil¬ 
icate” glass using a gemological refractometer is 
1.60±0.01, which is close to values for other glasses 
with like density. For the ultrathin slabs, inter¬ 
ference fringes occur near the Christiansen mini- 
mumj3 below, we use reflectivity data to compute 
n and ultrathin slab thickness. For the thin films, 
because the fringes are in the transparent region 
when n =1.6, we instead recorded the spacing of 
the diamond faces to estimate the film thickness. 
Spectra were aquired for an ultrathin slab with 
L =38/rm, two thick slabs with L =0.47 mm and 
2.364 mm, and two micron-scale thin films. Using 
n =1.6 gave 7.8/rm for the thick film and 3.5/im 
for the thin film. Later spectral scaling suggests 
that L =1.78 /im for the thinner film. The thicker 
of the two thin films was double the absorbance 
of the thin film, which is also consistent with the 
fringes. 


4.2. Laboratory spectroscopy 


We used several laboratory detectors at Wash¬ 
ington University in St. Louis to measure room 
temperature (18-19°C) spectra across the UV-IR 
wavelength range. In the UV-VIS (A =190-1100 
nm, A =0.19-1.1/tm, or v =52,630-9090 cm“^), 
we acquired unpolarized transmission and specu¬ 
lar reflectivity spectra using the double-beam UV- 
1800 Shimadzu@UV-VIS Spectrophotometer to 
probe the mid-UV and the UV tail at high res¬ 
olution (1 nm). The UV-VIS spectra were not 
acquired in an evacuated chamber. 


The infrared spectrum for “cos mic silicate” 
glass was previously published in 
(20H). T he same methodolo gy from 


S2ec^.et^ 


Speck et al 


( 2011 ) and Hofmeist^ ( 2003h was applied to mea¬ 
sure the IR spectrum for metallic iron in this work. 


^The Christiansen minimum occurs where the refractive in- 
dex of a medium changes rapidly and may approach that 
of the surrounding medium (usually air or v acuum). Th is 
results in very little reflection or absorption jCone 11 1 19691) . 
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4-.2.1. Spectral merging 

For the “cosmic silicate” glass, spectral seg¬ 
ments acquired for the two thin films, the ultrathin 
slab, and the two thick slabs were scaled in the 
region of wavelength overlap and merged. When 
scaling, the segment with the lowest absorbance 
was presumed correct. The lowest absorption that 
we can measure is the “noise” at the low point in A 
near A = 2.5 /im (4000 cm“^). We therefore shifted 
the spectrum up to avoid negative values; the 
smallest A is 0.001 or 0.0001, consistent with the 
signal-to-noise ratio of the spectral segment. We 
scaled the thin film spectra to match the upturn 
at low frequency (high wavelength) for the slab 
sample. We used the data over the Si-0 stretch¬ 
ing peak for the thin film, but otherwise used the 
thick film. Spectra from the thick (L ~0.47 mm) 
slab, rather than the ultrathin slab {L ^38^m), 
were used in the mid-IR due to noise considera¬ 
tions, although the data on the ultrathin slab were 
used to guide merging. 

For the reflectance spectra, scaling was not nec¬ 
essary; spectra were merged so that the noise was 
minimized. For Fe, reflectivity is nearly 100% in 
the far-IR; attainment of this limiting value in¬ 
dicates that the absolute values were recorded. 
For the “cosmic silicate” glass, the reflectivity in 
the visible is in accord with independent measure¬ 
ments of n. 


(damped harmonic oscillator model) assumes that 
the spectral peak breadths are due to damping 
(mode interactions), but in glasses, peak breadths 
also arise from a continuum of bond lengths in the 
glass and multiple overlapping peaks in the spec¬ 
trum. Likewise, Kramers-Kronig analysis fails to 
produce physically reason able values when peaks 
are u nresolved doublets (jOiesting fc Hofmeister 


200211 . 


We use an alternative, more direct method to 
derive n and k values for our “cosmic silicate” 
glass. Using absorption, reflectivity, and sample 
thickness, we determine the imaginary index of re¬ 
fraction k{X) via Eqs. [2]-[4]and then use our reflec¬ 
tivity spectra R to calculate n(A) using Eqs. HHS] 
At all A, k was determined from absorption data 
where we corrected for reflections. Let = 

transmittance, such that lab absorbance, in com¬ 
mon logs is represented by Eq. [51 


aSCM = —logio{Itrans/Io). 


( 2 ) 


For the polished slabs, we assumed that the 
measurement includes back reflections (R) from 
the slab’s top and bottom faces and compute 
the true or ideal absor ption coefficient ( A) after 
Hofmeister et al.l ( 2009ll and Speck ( 2014 


A = 2.3026[ascM + 2logio(l - R)l/L. (3) 


4.3. Optical function derivation 


We derived the optical functions for both “cos¬ 
mic silicate” and metallic iron which are shown in 
Figures 0] and [S] 


For the metallic iron spectra, we performed 
Kramers-Kronig analysis on the refl ectivity data 
to obt ai n n and k. For details, see e.g . . ISnitzer et al.l 
( 1962il : lAndermann. Caron fc Dons (Il965[) How- 
ever, applying the standard methods for our “cos¬ 
mic silicate” was problematic. Reflectivity R for 
Mg-Ca-Si-0 glasses is low (20% at the peaks, com¬ 
pared to ~80% for c rystalline silicates with siin - 
ilar composition; c.f. iMerzbacher fc Whitel 1988). 
Consequently, the standard methods used to de¬ 
rive n and k for crystalline silicates sometimes 
fail to provide positive values for the dielectric 
function €2 at all A for amorphous or glassy sili¬ 
cates. This problem can be understood in terms of 
glass structures. Specifically, classical dispersion 


For the “cosmic silicate” glass, Eq. [3] overcor¬ 
rects because the reflection from the back surface 
of the sample is less than R due to light being at¬ 
tenuated slightly. Thus, we adjusted the spectra 
manually for this effect and set i? = 0 in Eq. [31 for 
the “cosmic silicate” glass. We then calculated the 
imaginary part of the complex refractive index k 
using Eq. 01 


k\ = 


A 

47riy 


(4) 


The value of k is fixed at 0 for =0 cm ^ and 
for very large v above the X-ray region. 

Laboratory reflectivity (R) equals Iron- 

free glasses are very transparent and R is thus very 
low in the near-IR to visible; in this range, we 
computed R from Fresnel’s relationship, Eq. O 


R=[(n-l)^ + k^][{n+lf + k% (5) 
































using measured n„is and k =0. At other wave¬ 
lengths, we invert Eq.[5]to provide n(A) from mea¬ 
sured R and k (Eq. El) : 


1-R 


n{X) = z -5 + +(l-5) (6) 


1-R' 


For metallic iron, our n and k values extend to 
at A =0.19 ^m; we did not extrapolate to shorter 
wavelengths. All metals have nearly 100% reflec¬ 
tivity in the far-IR, w hich ar i ses fro m the nearly 
fr ee electron mod el of Drudd ( IQOOll as discussed 


WootenI (1197211 . 


For the “cosmic silicate” glass, because the 
composition was nearly 1 MgO -I- lSi02, n and 
k were averaged fro m data on MgO and silica 
glass in IPalik ( 1985ll . Also, n is pinned at unity 
above the X-ray region and k goes to 0 at very 
short wavelength. At long wavelengths (far-IR), 
we assumed that k decreases to 10“^ at 0.001 cm“^ 
(10 m and that n asymptotes to 3.16 at 0.001 cm“^ 
The value of k is an estimate, based on the value 
of k approaching 0 as the frequency goes to zero. 
The value of n is from the low frequency dielectric 
function, calculated below. 

According to Apoen fc Breskei ( 1952h . the di¬ 
electric constant is additive and can be estimated 
from the oxides using: 


e = Ecip* 



Wavelength (/.im) 


Fig. 3.— Comparison of complex refractive in¬ 
dices (n & k) for “Cosmic” and “Astronomical” 
silicates. Solid line = this work; dotted and dashed 
lines = past works. Insets are close-up views of 
the UV-vis region; Top: real part of complex re¬ 
fractive index; Bottom: imaginary part of complex 
refractive index; x-axis is wavelength in /rm. 


where pi is the fraction of the i*^ component and 
Ci was o btained at 0.45 GHz (0.01 5 cm“^) from the 


table in lAppen fc Breskeil (11952 1. See also the ta¬ 


ble on p. 318 of IScholze ( 199lH^ 


5. Comparing optical functions 

Our new optical functions for “cosmic sili¬ 
cate” glass and metallic iron are compared to 
popula r optic al functions cur rently in use, i.e. , 
Drainel ( 2003l l for silicate and Ordal et al. ( 1988h 
for metallic iron in Figures |4] and El Figure |4] 
clearly demonstrates that at A > 0.2 p,m the dif¬ 
ferences are huge. In particular, the imaginary 
part of the complex refractive index (fc) is several 


*also available here: http://www.lehigh.edu/imi/docs_GP/ 
Slides/GlassProp_Lecture22_Jainl.pdf page 24 


orders of magnitude lower in our new “cosmic sili- 
cate” g l ass me asurements than in those derived by 
Draine ( 2003l l. In part this difference comes from 


the difference in methods. The new data are deter- 
mi ned from exper imental laboratory data whereas 
the Draine ( 2003f) data are interpolated in this re¬ 
gion and were specifically fixed to get high enough 
opacity to match observations. 

Figure El shows that our new optical functions 
for metallic iron matc h the overall streng th and 
shape of the data from lOrdal et al.l ( 19881) . How¬ 


ever, our higher resolution spectroscopy measure¬ 
ments reveal details with the optical properties of 
metallic iron not seen in the Ordal data. In partic¬ 
ular the real part of the refractive index (n) drops 
to a very low value at A ^20 /xm and ~ 50 /xm; 
meanwhile, the imaginary part, k, shows a bump 
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Fig. 4.— New optical functions of “cosmic sili¬ 
cate” glass from this work (blac k lines), compared 
to “a stronomical silicate” from IPrainel (red lines 
200 , 11 . 



1 10 100 


Wavelength (//m) 


Fig. 5.— New optical functions of metallic iron 
from this work (black l ines), compared to values 
from Ordal et aD ( 1988ll (red lines). 


around 30 In addition, the new optical func¬ 
tions extend down into the UV-Visible part of the 
spectrum, rather than ending at the NIR/visible 
boundary. This is very important for modeling in 
the interstellar medium or stars above about K 
type, because the majority of the flux is emitted 
in this higher energy region. 


Our values differ from and bet ter represent Fe 
than those of lOrdal et al.l (|l988[ l for the follow¬ 


i ng re asons. The tables of n and k in lOrdal et al 
(jl988[ l are widely spaced (by 20cm“^) in the far- 
IR; the spacing is even greater at sh orter wave¬ 
lengths. Also, the Ordal etaP ( 1988ll data were 
obtained fro m a Kramers- Kronig analysis that in¬ 
corporated a lPrudd ( 190Clll model at high and low 
frequencies (short and long wavelengths). We did 
not use a Prude model. Kramers-Kronig analysis 
is effectively a smoothing function; this smooth- 
ing, combined with their low resolution, makes the 
lOrdal et al.l 1 19881) n and k data very smooth and 
prohibits resolution of narrow features. Our high- 
resolution iron data resolve the narrow features. 
In the visible, we also see a feature at ~0.4-0.5 
in our n and k d ata that wa s observed in previous 
studies (see e.g. |PaliPll99lll with a discontinuity 
between their non-overlapping data sets. Our data 
are continuous across this feature and so improve 
upon these past works. In the near-IR, the small 
feature at 2-3 fim in our iron data is not an arti¬ 
fact and is observed due to our high instrumental 
resolution. Therefore, we have provided improved 
n and k at higher resolution. The trade off is the 
noise in the far-IR longer than 10 /rm; however, as 
will be seen in §|6l this noise does not affect our 
radiative transfer models of the observations. 


6. Radiative transfer modeling 


6.1. DUSTY 

We use d the 1-P radiative t r ansfer program 


DUSTY (Ivezic fc Elitzur 1995: Nenkova et al 


2 OOC 1 II to determine how our new “cosmic sili¬ 
cate” glass optical functions perform against pre¬ 
vious “astronomical” and “circumstellar” silicates 
n and k. For our test case, we built simple models 
for HP 161796, a dusty oxygen-rich pre-planetary 
nebula/post asymptotic giant branch star. We 
chose HP 161796 for several reasons: (1) it has 
been observed many times over several decades 
and over a large wavelength range, and thus has a 
well sampled spectral energy distribution (SEP), 
as shown in Figure [6] (2) its central star is hot 
enough to emit significant UV radiation, but not 
so hot that photo-dissociation and ionization play 
a major role in the nebula, simplifying the mod¬ 
eling; and (3) although not co mpletely spherical , 
it is a remarkably round object ( Ueta et al. 1999ll 
and thus using a 1-d RT code is a reasonable ap¬ 
proximation. 
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Fig. 6.— Spectral Energy Distribution (SED) of HD 161796. Very small black dots are the ISO SWS data; 
very small grey dots are the ISO LWS data; medium black filled circles are the short wavelength IRAS LRS 
data; medium grey filled circles are the long wavelength IRAS LRS data; large filled circles are the IRAS 
photometry data (grey are uncorrected; black are color corrected). Sources for other photometry points are 
given in the legend. a;-axis is wavelength in microns; y-axis is flux (AE\) in W 


6.2. Modeling HD 161796 with synthetic 
optical fnnctions 

We have generated models using the new “cos¬ 
mic silicate” glass n, k data presented herein as 
well as using synt h etic o ptical f unctio ns from 
Draind (l2003l) : |di] (Il984l) : IOHM] (Il992ll . which 
are compared in Fig. [T] 

For our models we adopt the central star tem¬ 
peratu re, T* = 6750 K, from iHoozgaad et al 
(120021) . which falls within the currently accept¬ 


able range for the object’s central F3I b class star 
(6699 K. 6600 K. 6700 K. 6900-7350K: iKovtvukb 
2007; Lvubimkov et al.lI 2 OIOI : Land 19921 respec ¬ 


tively) and previous models (ISkinner et al.l 11994 : 
Meixner et al. 20021) . Fig. [7] shows that the effect 


of v arying this A valu e by ±150 K is negligible 
(c.f. DePew et al. 20061) . DUSTY and many other 
RT models invoke Mie theory to calculate absorp¬ 
tion and scattering cross-sections when a user sup¬ 
plies n and k. The simple models shown here are 
made with the default DUSTY settings, which im¬ 
plicitly assume separate populations of spherical 
grains. However, spherical grains are probably in¬ 
appropriate in many astrophysical environments 
and are not the best representation of silicate 
features or of metals, whose FIR behavior is ex¬ 
tremely shape and agglomeration dependent. For 
’’real” models, we recommend supplying absorp¬ 
tion and scattering coefficients for non-spherical 
grains calculated directly from n and k. In that 
way, our new n and k values allow users to in- 
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vestigate the effect of grain size an d shape (e.g. 
Pollack et al .11 1993 : Min et al. l2003l l. 


Figure [7] (top) sho ws the effect o f using the best- 
fitting model for the Draind ( 2003li function m odel 
with t he original “astronomical silicate” from IDL 
and with the warm and cold “circumstel- 


lar silicate” optical functions from lOHMl (Il992ll . 
Figure [7] {bottom) shows the best-fitting models 
using the origi nal “ astron omical silicate” optical 
functions from |DL| ( 1984 ) and those from lOHM 
(Il992r) . all three of which are embedded options 
in DUSTY. All model parameters are shown in 
Tab l e El The difference in p erformance between 
[Pll (jl^4h and IPraind (|2003 ll is that the V-band 
optical de pth is very sl ightly lower when using 


Draind ( 2003 1 optical function. Using 
1992h synthetic functions requires even 


the newer 
the lOHMl 
lower V-band optical depths. 

We assumed a radial dust density distribution 
of to reflect a constant mass-loss rate. We 
also initially assumed that the outer radius of the 
dust shell (i?out) is 100 x the inner radius (i?in), 
although the models are not ver y senesitive to 
this parameter. ISoeck et~ ( 20091 and references 
therein) demonstrated that the inherent degen¬ 
eracy in RT modeling means that this number 
can be changed by compensating using the inner 
dust radius/temperature. Therefore we arbitrar¬ 
ily chose this value for comparison purposes. We 
used an MRN grain-size distri bution to be consis¬ 
tent w i th the assum p tions of Draind (12003 1: DL 
( 1984 1: OHM ( 1992 1. Speck et al.l ( 20(1911 showed 
that DUSTY models of spectral energy distribu¬ 
tions are not very sensitive to the precise grains 
size distribution as long as the distribution does 
not exclude many small grains or include very big 
grains (>1 /im). 


Table 2: Best-fitting DUS'T V models u si ng syn¬ 
thetic optical funct ions from IDLI (|l984h : iDraine 
(l2003ll:[OHh3 (Il992ll 


source 

TV 


D2003 

2.25 

125 

DL84 

2.35 

125 

OHM92-C 

1.75 

125 

OHM92-W 

1.75 

115 


6.3. Modeling HD 161796 with our new 
“cosmic silicate” optical functions 


Our purpose in developing these simple models 
(an example of which is shown in Fig. [5]) is not to 
make an exhaustive study of the parameter space, 
or to optimize the fit to this particular object’s 
SED, but instead to demonstrate the effect of us¬ 
ing the new, real mineral data presented herein 
compared to more commonly use “synthetic” com¬ 
plex refractive indices. We have deliberately cho¬ 
sen the minimum number of dust components nec¬ 
essary for a fit and commonly used model param¬ 
eters such as the ^ radial density distribution 
and the MRN grain-size distribution with spheri¬ 
cal grains. The effect on radiative transfer mod- 
eling of choosing o ther parameters is discussed by 


Speck et al. ( 20091 1. 


6.3.1. “Cosmic silicate” ONLY 


The initial models assumed a single dust com¬ 
ponent; later models included metallic iron and 
“cosmic silicate”. We also included external heat¬ 
ing of the dust shell by the interstellar radia¬ 
tion field as well as a pile-up of material in the 
outer regions of the dust shell where the stellar 
outflow meets the interstellar medium in order 
to investigate whether composition or other fac- 


tors 

improve the match to the shape of the SED 

(c.f. 

Gillett et al. 19861: Young. Phillips &: Knannl 

1993 

). 


Models using only our new “cosmic silicate” 
glass optical functions do not match the data 
or the models using synthetic optical functions. 
There are two issues: (1) The shape of the SED 
produced using “cosmic silicate” is too symmetric 
and cannot match the observed emission at far- 
IR wavelengths; and (2) the visual optical depth 
required is unfeasibly high (100 < Ty < 150). 
We investigated the effect of varying several model 
parameters in an attempt to alleviate the far-lR 
problem including the radial density distribution, 
the geometrical size of the dust shell, inclusion of 
crystalline silicates, inclusion of a step function 
in the density distribution to replicate the poten¬ 
tial pile-up where the circumstellar materials meet 
the ISM, and inclusion of external heating of the 
dust shell by the inter s tellar radiation field (see. 


Gillett et al. 19861 : Young. Phillips fc Knannl 


19931). None of these resolved either problem with 
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Fig. 7.— Best-fitting DUSTY models for “Astro¬ 
nomical Silicate”. a;-axis is wavelength in microns; 
y-axis is flux (AFa) in Wm“^. Main panel is a 
linear-linear plot; inset is a log-log plot of the same 
data. In all panels, large black filled circles = pho¬ 
tometry points from Figure 1; very small dots = 
ISO spectrum. Top: Best model fit is achieved us - 
ing complex refractive index from Draind ( 2003ll . 
The effect of varying T* by ±150K is negligible. 
Bottom: Best-fitting DUS TY models using co m- 
plex refractive index from IPraine ; [dI]: IOHMI for 
Tatar =6750K and varying opacities. 


the “cosmic silicate” only models. 

6.3.2. “Cosmic silicate” and metallic iron mod¬ 
els 

As discussed above (§ HD, grains containing a 
mixture of metallic iron and iron-free silicates have 
long been invoked to account for the low opacity 


of sil i cates either explicitly (e.g. iMerrill fc Stein 
1976 : OHMl[l992h or implicitly (e.g. Draine 200,11 
and references therein). Condensation models sug¬ 
gest that iron should form as metallic iron grains 
rather than as part of a solid solution of ferro - 
magnesian silicates (e.g. Fodders fc Feelev 1999ll . 
Therefore we also included metallic iron grains 
into our models, along with “cosmic silicate” glass, 
in varying proportions. DUSTY assumes that the 
grains are co-spatial; i.e. the absorption and scat¬ 
tering cross-sections are calculated assuming the 
different grain compositions have the same size 
distribution and have the same temperature (es¬ 
sentially it a composite grain, comprising a silicate 
sphere and an iron sphere in thermal contact.) 

Figure E] shows that a modest addition of metal¬ 
lic iron dramatically improves both the shape and 
the optical depth requirements for the models. 
With only 25% metallic iron grains, we can re¬ 
duce the optical depth to ^2, comparable to the 
models with synthetic optical functions. 

Whereas including iron into the models im¬ 
proves the values for the visual optical depth, we 
still have a modeled SED that is narrower than 
either the observations or the models using syn¬ 
thetic optical functions. Inclusion of crystalline 
silicates (which are obviously present based on the 
ISO spectrum in Fig. E] and have been includ ed in 
previous models, e.g. Hoozgaad et al.l 2002li also 
does not alleviate the problem of SED shape and 
produces short wavelength crystalline features not 
seen in the observational data. 


6.3.3. Two-layer models 

Consequently, we used DUSTY to build a 2- 
layered model, which is shown schematically in 
Fig. O The inner part of the shell is modeled 
to get the best fit shortward of 35/xm. The 
output from this model is then the inner cen¬ 
tral source for the second layer model. Because 
RT modeling is inherently degenerate, there are 
several models that produce good fits, summa¬ 
rized in Eig [9l All good fits consist of an in¬ 
ner region dominated by “cosmic silicate” glass, 
with 25% metallic iron, and a geometrically thin 
outer zone where the composition no longer con¬ 
tains any glassy silicate. Instead the thin outer 
layer contai ns 10% metallic iron, 30% Fe-ri ch py - 
roxene from iHenning fc Mutschke (bronzit e; 1997) 
and 60% forsterite ffrom IPitman et al. 2013aD) . 
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Fig. 8.— Best-fitting DUSTY models for “cos¬ 
mic silicate” with metallic iron included, cc-axis 
is wavelength in microns; y-axis is flux {XFx) in 
Wm“^. Main panel is a linear-linear plot; inset 
is a log-log plot of the same data. In both pan¬ 
els, large black filled circles = photometry points 
from Figure 1; very small dots = ISO spectrum. 
Thin black dashed line = model wi th “cosmic sil- 
icate” n &L k from this work, but iDraind (12003^ 
best-fitting model paramete rs; thin g rey so lid line 
= best-fitting model using lOraind ( 2003 ): thick 
black solid line = best-fitting using “cosmic sili¬ 
cate” n & /c with 10% metallic iron; thick black 
dotted line = best-fitting using “cosmic silicate” 
n Sz k with 20% metallic iron. 


Henning fc Mutschkd (1997) do not state whether 


the iron (FeO) content is given as mole % or 
weight%. However, in either case the sample 
should be referred to as Fe-rich enstatite rather 
than bronzite. We use this mineral as an example 
of a single crystalline silicate, rather than an at¬ 
tempt to match all the features. The best 2 layer 
model is shown in Fig. [TUI 


This model may suggest that crystallization of 
silicates occurs either as a result of shock heating 
of grains where the circumstellar medium ploughs 
into the surrounding ISM; or as a result of high 
energy (FUV) photons coming from the ISRF. It 
is possible that this thin layer is also where many 
of the silicate grains are destroyed, which would 
account for the difference between circumstellar 
and interstellar silicate features. Indeed, it has 
been posited that most interstellar silicates form 


in cold, dark molecular clouds rather than being 
ejected from dying stars (e.g., IZhukovskal 12008 : 


Drainell2009l) . However, this idea is controversial 


(see, e.g. [Jones fc Nuthll201ll) . 



Fig. 9.— Schematic cartoon of the 2-layer DUSTY 
models 


6.4. Iron optical functions 

There is a subtle effect as a result of replac¬ 
ing the widely used metall ic iron optical func¬ 
tions from Ordal et al. ( 19881) with our new metal¬ 
lic iron optical f unctions. In order to match the 
model using the Ordal et ahl ( 19881) data we had 
to either increase the visual optical depth or in¬ 
crease the fractional abundance of metallic iron in 
the model. 


7. Discussion 

It is clear from the way in which the optical 
functions for “astronomical silicates” were derived 
that they have some limitations. At the short 
wavelength end, the materials used do not actually 
represent the solids we expect to be dominating as- 
trophysical environments. At the long-wavelength 
end, the data are derived from astronomical ob¬ 
servations and assumptions and thus provide no 
flexibility to match environments with different 
dust components. Past attempts to derive the op¬ 
tical functions for astrophysically-relevant silicates 
were hampered by the lack of availablity of consis¬ 
tent laboratory data at that time. Decades later. 
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Fig. 10.— Best-fitting DUSTY models using two- 
layer approach, cc-axis is wavelength in microns; 
y-axis is flux {XF\) in Wm“^. Main panel is a 
linear-linear plot; inset is a log-log plot of the same 
data. In both panels, large black filled circles = 
photometry points from Figure [51 very small dots 
= ISO spectrum. Thin black dashed l i ne = model 
with cosmic silicate n Si k, and lPraine ( 2003 1 best¬ 
fitting model parameter s; thin g rey so lid line = 
best-fitting model using Draind ( 2003 1: thin red 
dotted line = inner shell model thick blue solid line 
= 2-layer model = inner model with thin shell of 
crystalline silicate and metallic iron on the outside. 
See schematic cartoon in Fig El 


we now have the means to overcome at least some 
of these issues. 


The 10 /im silicate feature occurs in many 
dusty astrophysical environments. It varies 
from object to object and even within a sin- 
gle o b ject both tempora lly 




Monnier et al.l 


19981 : ICuha Nivoei et ah 2011a) and spatially 
(e.g., in n Car, N. Smith, Pers. Comm; and 
ICulia Nivogi. Speck fc Volkiboilb l. Within a sin¬ 
gle type of astrophysical object, the feature shows 
huge variations in terms of its peak position, width 
and its ratio to t he ^ 18^m feature (e.g., ISoeck 
19981 : OHMI [1993) . Variations in feature or SED 


shape from object to object cannot be analyzed 
modeling with synthetic optical functions. 


Small circumstellar dust particles are likely 
much more compositionally and structurally in¬ 
homogeneous than the simple glass + Fe cases 
we have considered here. Thus, we are provid¬ 
ing data on well-characterized endmember glass 


to allow modelers to mix different compositions 
and impurities as they choose. Our data also can 
aid in distinguishing between dust that is glassy, 
crystalline, or some combination of the two, which 
is important because the degree of structural ho¬ 
mogeneity has implications for its formation, and 
subsequent proces sing, evolution and destruction 
( Speck et alJ 2011 1. 


7.1. Comparison to previous models 

HD 161796 is a well-studied object and many 
previous m odels exist, whi c h mak es it an excellent 
test case. Hrivnak et al.l ( 19891) produced mod¬ 
els of several candidate pre-planetary nebulae, in- 
cludi ng HD 161796 , using dust opacities derived 
from IVolk fc KwokI ( 19881) . Their output param¬ 
eters depend on distance, but are consistent with 
the later models in terms of inner dust radius and 
stellar luminosity. 


Justtanont et al.l (Il992i l made a very simple 


model of HD 161796’s SED using a 6000 L© F- 
type star surround by dust with a single blackbody 
at 90 K modified by the emissivity from disor¬ 
dered oliyine from iKratschmer fc HuffmanI (197^ 
and assuming an MRN size distribution. Their 
best-fitting model did not match the SED at all 
wayelengths, and although it was somewhat im¬ 
proved by the addition of magnetite, there re¬ 
mained mismatches between the model and the 
observed SED. 

Skinner et al.l (119941) used multiwavelength 


imaging to constrain their models of HD 161796. 
Their model assumed an effective star tempera- 
ture of 630 0 K and used astronomical silicate from 
|DL| (1984) in an MRN distribution. Their best¬ 
fitting model gave a stellar luminosity of 36OOL0 
and an inner dust radius of 8.9x10^^ cm. 

Up to this point, the past models assumed 
spherical symmetry (as do we). Furthermore, 
there was not yet observational evidence for crys¬ 
talline silicates. With the advent of ISO, more 
detail in the far-IR spectrum was revealed and 
led to modeling incl uding crystalline silicates. 
Hoozgaad et al. ( 20021) modelled the full 2-200 ^m 
ISO SWS -I- LWS spectrum together with addi¬ 
tional data extending the SED from the optical 
to the millimeter wavelength range. Their model 
assumed a star of 3OOOL0 at 6750 K with an inner 
dust radius of 21x10^® cm (and outer radius of 
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63x10^® cm). iHoozgaad et al. ( 20021) developed a 
two layer model in which the inner dust is dom¬ 
inated by amorphous silicate and the outer dust 
includes a mantle of water ice. They used a num¬ 
ber of different optical functions in their model in 
order to cover the full wavelength r ange, but the 
amorphous silicate was solely from lOHMl (Il992h 
and was extrapolated out to the extremes of the 
wavelength range. In addition, they included crys¬ 
talline water ice, as well as crystalline silicates in 
the form of enstatite and forsterite. They fo¬ 
cussed on developing two parts of their model - 
the core-mantle grains and the crystalline silicate 
fraction, although they found that the lack of sen¬ 
sitivity of the models to the precise abundances of 
crystalline silicates made abundance constraints 
weak, but found that ^^10% by mass was reason¬ 
able. The average temperature for the amorphous 
silicate grains was found to be 110 K (and ranged 
from 47-137 K) while the crystalline grains were 
cooler (at ~50-90K), consistent with our results. 

Meixner et al.l ( 20021) produced a 2-d model of 


HD 161796 for which they combined both mul¬ 
tiwavelength imaging and the SED to constrain 
their models. Their best-fitting model had a 
2800 Lq star with an effective temperature of 
7000 K. They found an inner dust temperature 
of 110 K and a corresponding inner dust radius of 
11x10^® cm. Whereas the model is axisymmetric, 
the change in dust density from equator to pole 
is small. Their model included both forsterite 
and enstatite but the coarse spectral resolution 
precluded fitting the narrow crystalline silicate 
features in detail. The grain size distribution 
is KMH-like with a higher upper limit on grain 
size, and they assumed the dust composition from 


Hoozgaad et al. ( 20021) . 


Given the range of optical functions and model 
types used, it is interesting (and reassuring) that 
all models including ours end up being dominated 
by amorphous silicate dust at ~100K. However, 
previous models used synthetic optical functions 
for amorphous silicates and thus did not include 
the correct Fe/Mg fraction and cannot reveal the 
shell structure seen in our modeling results. 


8. Conclusions 

We have presented new optical functions for 
iron-free glassy silicate with chondritic compo¬ 


sition and metallic iron, extending from UV to 
far-IR wavelengths (0.19-250/rm), which are mea¬ 
sured directly in the laboratory using the same 
sample throughout. The “cosmic silicate” data 
extend to the X-ray region. These new opti¬ 
cal functions improve upon prior synthetic opti¬ 
cal functions for “astronomical” and “circumstel- 
lar” silicates by preserving the chemical composi¬ 
tion, sample properties and the method of optical 
function derivation (all lab data). A simple model 
excursion using HD 161796 as an example shows 
that: 

1 . using synthetic optical functions automati¬ 
cally fits most SED slopes and thus provides 
false parameters for models. 

2 . inserting back-calculated n and k from ob¬ 
served data is not necessary to achieve the 
desired SED slope; only a small amount of 
metallic iron is needed to make the models 
sufficiently opaque; 

3. the precise ratio of iron/cosmic silicate 
strongly affects the far-IR emissivity law. 
Thus, our new laboratory-derived opti¬ 
cal functions can be used to determine 
how composition affects observed long- 
wavelength emissivities and emissivity in¬ 
dices (/3 value), which is relevant to under¬ 
standing many long wavelength observations 
including those from ALMA and Herschel, 
Spitzer, IRAS and ISO] 

4. our models of HD 161796 hints at crystal¬ 
lization of silicates at the outer edge of the 
circumstellar shells. While this is a prelimi¬ 
nary result that requires further study, this 
suggestion is not apparent when using syn¬ 
thetic optical functions; 

Because our cosmic silicate has chondritic 
abundance, it may not be precisely the composi¬ 
tion of glassy silicates formed away from the solar 
neighborhood. Future work will include other sil¬ 
icate glass compositions including both iron-free 
and iron-bearing compositions. 

This material is based upon work supported 
by the National Science Foundation under Grant 
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